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ABSTRACT 

Via the Rossiter-McLaughlin effect, it is possible to measure the sky-projected angle between the stellar spin and a planet's orbital 
spin. Observed orbital inclinations have been found to range over all possible angles. A tentative detection of a correlation between 
the dispersion in spin/orbit angle and the youth of the system is revealed, using spin/orbit measurements for hot Jupiters around stars 
with masses > 1.2 Mq for which age estimates are more accurately determined. The chance of this pattern arising by chance has been 
computed to 7%. This appears in accordance with tidal dissipation where non-coplanar hot Jupiters' orbits tidally realign. The results 
show they would do so within about 2.5 Gyr. For the considered sample, the results give support to hot Jupiters being placed on non 
coplanar orbits early in their history rather than this happening late. Such events could involve strong planet-planet scattering. 
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1. Introduction 

For transiting planets, the Rossiter-McLaughlin effect (|Holt| 



Matsumura et al. 



1893^ Rossiter 



Gaudi & Winn 



[T9241 |McLaugh Un|fT924} |Queloz et al.|| 2000[ 
2007| l, allows the measure of /3 (also called A in 



the literature), which is the projection on the sky of the obliquity 
iff, between the stellar spin axis and the orbital spin axis. 

Up until recently planets were thought to be mostly on orbits 
coplanar with their star's equator (F abrycky & Winn 2009) 



something in line with predictions of disc migration ( |Lin et al. 



1996 Ward 1997 1. More recently a number of papers have 



shown that hot Jupiters on non coplanar orbits are common. 



2008, ^Moutou et al. 



including some planets on retrograde orbits (Hebrard et al. 

20091 [Narita et al.|. 2009b. Winn et al. 
I2OI0I IQueloz et al.||201o( [Triaud et al. 



2009a) [Anderson et a 



2010 1. Those measurements have been interpreted as showing 
that dynamical events are probably not uncommon and that not 
all systems can be understood by disc migration alone. Strong 



dynamical events such as planet-planet scattering (Rasio & 


Ford[[1996 


[Juric & Tremaine[[2008 Chatterjee et al.[[2008 1, or 


more secu 


ar processes such as Kozai-Lidov oscillations (Wu[ 


et al.|2007 


Fabrycky & Tremaine|2007| Nagasawa et al. 


2008 


Naoz et al. 


201 1)1, or chaotic interactions (|Wu & Lithwick 


201 If 



would place a planet on a highly eccentric orbit, whose passage 
at periastron is sufficiently close that tidal dissipation causes the 
planet to lose angular momentum and circularise around its star 

Understanding the origin of hot Jupiters is one of the keys 
to shedding light onto the processes that act during planet 
formation as well as those acting after planets have formed. 
Those processes allow us to place constraints on what happened 
and did not happen in our own Solar System. They will also 
help us match more accurately theoretical predictions of planet 
formation done in population synthesis simulations to the 
parameter space that planets currently occupy, as given by the 



observations (eg. |Ida & Lin| ( |2004| and [Mordasini et al.[p009[ )) 
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( |2010 a^ remark that if misaligned hot 
Jupiters do not require disc migration, aligned planets are not in 
contradiction with a scenario involving dynamical interactions 
and tidal migration, as planets will tend to realign with the star 
(see also [Hut| ( [T98T] l and [Barker & Ogilvie[ ( |2009l l). 

'Winn et al. (|2010a) point out a correlation between the 
stellar effective temperature and the spin/orbit angle. For stars 
with Teff > 6250 K, fewer aligned systems are found compared 
to stars with lower effective temperatures. This would show that 
tidal realignment timescales are different for different stars, as 
proposed by [ZaiEi| ( |i977) l in the context of binaries. Schlaufmanj 
( [2010 1) presents an independent confirmation of that correlation, 
using a different methodology. 

The aim of this letter is to combine the observational facts 
and offer an explanation. The results will then be discuss in light 
of the currently available theoretical framework. 



2. Motivation 

The lack of aligned systems for stars with T^ff > 6250 K that is 
noticed in [Winn et aT] ( |2010a) could also be explained by stel- 
lar physics combined with an observational bias: as predicted by 
stellar evolution, stars with masses greater than about 1 .2 Mq 
start on the Zero Age Main Sequence with temperatures higher 
than 6250 K. When H-core burning stops, they have cooled by 
several hundred Kelvin (fig.[T]). They do so in 3 to 4 Gyrs. This 
means that, while the planet and the star progressively realign, 
the star itself cools down. We are thus left with an aligned planet 
around an older, cooler star. Some, more massive, stars will cool 
to temperatures above 6250 K, but the timescale for realignment 
might be longer than the Main Sequence lifetime. Once they 
leave the Main Sequence, stars becomes too large for planets 
to be discovered by ground-based transit surveys as the contrast 
becomes too small. We thus have a bias to see misaligned plan- 
ets around hot stars, notably because we may not detect their 
aligned population. 
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Fig. 1. Main Sequence showing the Geneva stellar evolution 
tracks for solar metallicity as presented in Mowlavi et al. (sub- 
mitted) and plotted using (in Rq) as a function of T^-g. 
Tracks are labelled in units of Mq. Dashed line show the 2 
Gyr isochrone. Overplotted are the systems for which we have 
Rossiter-McLaughlin measurements. Aligned systems are red 
circles, misaligned systems are blue triangles. Higher metal- 
licities will move the tracks to the right. Data obtained from 

EXOPLANET.EU 



This explanation could be combined to the different realign- 
ment timescales described in |Winn et al.| ( |2010a( ) and |Zahn| 
1977| l since, as the star ages and cools, its convective zone 



would become larger too. If that explanation is right, we should 
expect a correlation between stellar age and alignment. 

The average stellar density, p*, is obtained directly from 
the planetary transit signal ( [Sozzetti et al.|[2007| l, the effec- 
tive temperature, Tf-g, and metallicity, Z, can be obtained via 
spectral analysis. Stellar mass and stellar age can be estimated 
from interpolating the stellar evolution tracks in (p*,reff,Z) 
space. Interestingly, stars > 1 .2 Mq spend less time on the Main 
Sequence, but increase their radii more than solar mass stars do. 
We thus have a higher resolution on the tracks to estimate ages 
on more massive stars than on solar mass stars. Such a subsam- 
ple should give the most precise and accurate ages that we can 
get. This is the sample used in this letter 
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Fig. 2. Secure, absolute values of /3 against stellar age (in Gyr), 
for stars with M^. > 1.2 Mq. Size of the symbols scales with 
planet mass. In blue squares, stars with M* > 1.3 M©; in red di- 
amonds 1.3 > > 1.2 Mq. Horizontal dotted line show where 
aligned systems are. Vertical dotted line shows the age at which 
where misaligned planets start to disappear. 



3. Sample selection 

Let us take only the most secure measurements for the projected 
spin/orbit angl^ for planets with stars > 1.2 Mq. There are 22 
objects in the sample (table [T}. The sample is divided in two: 
stars > 1.3 Mq (8 stars) and stars between 1.2 and 1.3 Mq (14 
stars). The angle and age estimates were obtained from the liter- 
ature, but for WASP- 17, whose error bar on the age was large. It 
was re-estimated for this letter, using the stellar parameters and 



density presented in Triaud et al. 



(2010 1 and interpolating in the 



Geneva tracks (Mowlavi et al. submitted). The new age estimate 
is 2.3 ±0.6 Gyr Its error bar is consistent with age measurements 
made by other teams. The new value is presented along with all 
other values in table [T] 

Plotting the absolute values of the measured projected 
spin/orbit angle /3 against stellar age (fig. [2|, a pattern is obvi- 
ous and as sharp as that presented in Winn et al. (2010a). While 
observationally, there should be no bias to preferentially detect 
aligned systems instead of misaligned systems at any age, stars 
older than ~ 2.5 Gyr show mostly aligned systems (rms - 22°, 
median = 5°). For stars that are younger we have a large range 
of obliquities (rms = 66°, median - 60°). Figure |3] displays the 
cumulative distributions on either side of the 2.5 Gyr age limit. 

To test the robustness of the pattern, a Monte Carlo simu- 
lation was performed taking the data with ages < 2.5 Gyr as a 
fiducial zone from which random samples of 8 measurements 
were drawn, allowing for repetitions. There is < 4% chance to 
draw a sample with median < 10° and rms < 60° which would 
allow a sample having seven aligned systems and one retro- 
grade system. If restricting the rms within 30°, similar to that 



Some measurements have been omitted for the following reasons: 
CoRoT-3 (sampling is poor dTriaud et al. 2009|> ), CoRoT-1 1 and Kepler- 

I Jenkins etal.|2010^ ) and 



8 (transits are incomplete ( Gandolfi et al. 



zora 



WASP-1 (angle is unsure ( jAlbrecht et al.|201 



0. 



Triaud A. H. M. J.: Age and spin/orbit angle 



3 



50 



(°) 



Fig. 3. Cumulative distributions in orbital inclinations for sys- 
tems younger than 2.5 Gyr (dashed blue), and older (plain red). 
For comparison, a uniform distribution (dotted black). 



observed, there is a probability < 1% that the distributions on 
either side of the 2.5 Gyr age are the same. Drawing randomly 
from the overall sample, there is a 2.6% chance to obtain a 
cluster containing 7 aligned systems and another at any angle 
> 20°. In addition a Kolmogorov-Smimov test was carried out, 
also comparing the distribution in p on either side of the 2.5 Gyr 
limit. AD- 0.661 is obtained corresponding to a probability 
of 1.2% that both distributions are the samq3 The same test 
shows that the distribution of angles around stars younger than 
2.5 Gyr has about 22% chance to be compatible with a uniform 
distribution, while for the older sample, this chance is of order 
10"^. By rearranging the data, selecting various cut-off and 
computing the KS test at each step, the probability of having 
two such different populations arising by chance is estimated 
to about 7%. It can be affirmed there is tentative evidence of a 
pattern in the data. 

We see that stars with masses > 1.3 M© are all younger than 
3 Gyr. Thus, when observing few ahgned systems on stars 
with Teff > 6250 A", Winn et al. ( 2010a i were in fact detecting 
an effect due to stellar age, or rather, time since planet formation. 



Like for all multivariate problems, figure |2] offers an incom- 
plete picture: it only shows two quantities in relation with time. 
At the moment orbital separations and mass ratios are quite sim- 
ilar since the bulk of the discoveries have been done by ground- 
based transit search programs. With increasing numbers of mea- 
surements over a larger parameter space we will eventually need 
to account for those extra parameters. 

4. Discussion 

The large variety of angles around the younger stars suggests that 
some misaligning mechanism happens during the youth of plan- 
etary systems. Notably, in combination with results by (Watson 
|et al.| ( |201 l| l showing no evidence for misaligned protoplane- 
tary discs, it lends strong support to a planet-planet scattering 
scenario occurring during the last stages of planet formation 
or soon in the aftermath of the disc dispersal like described in 
[Matsumura et al.| ( |2010bl l. 

When preparing figure|2] reason dictated that a dearth of old, 
misaligned systems was expected, not an absence. The complete 



lack of misaligned planets orbiting stars older than 2.5 Gyr in 
the current sample came somewhat as a surprise as secular in- 
teractions could place planets on inclined orbits well after the 
disc dissipated. A system presenting such characteristics can 
be found among the "older" systems: HAT-P-13, whose cur- 
rent configuration may have originated from secular interactions 
( |Mardling|20 1 0[ ) . If that history is right, its observed coplanarity 
may be a chance alignment. Chance alignments can occur eas- 
ily since firstly, we observe a projected angle, and not the 
real obliquity ^ and secondly, theoretical predictions such as Wu 
et al.|(2007] ), Fabrycky & Tremaine (2007) and Nagasawa et al. 
(2008 ) predict very high orbital inclinations, but also a number 
of aligned systems. 

There is great interest in matching those theoretical distri- 
butions to observations (notably for young hot Jupiters), but the 
evolving nature of the spin/orbit angle distribution makes this a 
tricky task. Multi-body dynamics are less concerned about abso- 
lute masses than about mass ratio. In systems where no Jupiter 
has formed, we would expect planet-planet scattering between 
Neptune-mass planets producing an inclined hot Neptune popu- 
lation. If the inital stages will be similar, the later ones will not: 
tidal circularisation and realignment timescales will be different. 
Spin/orbit angles for planets of masses < 0.1 Mjup will be less 
affected by tidal realignment and offer a closer picture of the ini- 
tial spin/orbit angle distribution than hot Jupiters. A hot Neptune, 



Hat-P-1 1 b has been recently detected misaligned by, Winn et al. 
( 2010c| l and confirmed by Sanchis-Ojeda & Winn ( 201 I| l. 

This work has focused on stars with masses > 1.2 Mq. If age 
is what determines primarily whether a hot Jupiter is observed 
aligned or misaligned, since solar mass stars are detected in 
average older than more massive stars, it is not surprising that 
their planets are coplanar. There nevertheless is an interest in 
looking at that population carefully which stems from work by 
Burkert & Ida| ( pUOT) , |Currie| ( |2059] l and Alibert et al.| (|20TT) 
who argue that discs around the more massive stars are not 
long lived enough to produce an aligned hot Jupiter population 
via disc migration. In the mean time, if planet formation is 
more efficient in more massive discs (found around more 
massive stars), then one could expect a higher occurrence of 
planet-planet scattering around such stars. If this is true, it could 
point towards two pathways for bringing hot Jupiters to their 
observed location which would be dependent on stellar mass. 
Unfortunately stellar ages are less precisely determined for solar 
mass stars as illustrated by the isochrone on figure [T] 

The change in the shape of the distribution of spin/orbit an- 
gles with time is indicative of some orbital evolution, presum- 
ably through tidal interactions between the star and the planet. 
Barker & Ogilvie ( 2009) show that retrograde planets decay 
into their star on timescales two to three times shorter than pro- 
grade planets would do, for given initial conditions. Their infall 
timescale for a typical, retrograde, hot Jupiter are of order of 
a few Gyrs. Winn et al. ( 2010aj ) present similar behaviour In 
addition they show that, for a given stellar mass, a more mas- 
sive planet will realign and in-spiral faster than a lighter on^ In 
both papers the retrograde planets realign with the star but only 
shortly before falling into it. It would thus be unlikely to observe 
them at these very particular phases. Nevertheless such exam- 
ples could be found in WASP-12, 18 and 19 (eg. [Helher efaT 



the same test on the pattern presented in 



Winn et al. 



1 2010a I gives 



6.1% chance that the distributions on either side of 6250 K are the same 



^ Incidentally, this could explain a second observed feature in relation 
to the angle p. As shown in |Moutou et al.| ( (201 1) , there is a lack of ret- 
rograde massive planets (> 5 Mj^p), something expected if retrograde, 
massive planets realign faster to their star than other planets do. 
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( |201 l| l). |Matsumura et al.| ( [2010a| l describe how planets initially 
placed on mildly inclined or aligned orbits, are less likely to in- 
fall and more likely to survive until observed. Nevertheless, in 
most cases tidal realignment corresponds to the disappearance 
of the planet. 

If retrograde planets plunge into their star as they tidally 
realign, a decreasing number of hot Jupiters should be observed 
with time. No such decreasing trend can be found when con- 
sidering all the hot Jupiters presented in the literature around 
stars in the mass range considered for this paper This is at odds 
with evidence of a trend between semimajor axis and stellar 
age showing a lack of very short orbits around older stars as 
presented by Jackson et al. ( 2009 ) who interpreted it as evidence 
of the destructive tidal orbital decay of hot Jupiters. Looking 
at the semimajor axes of the targets in table [T| a similar trend 
appears. This may indicate we do not have enough objects yet 
to detect the expected decreasing fraction of hot Jupiters with 
time. Alternatively the results presented here, and those of |Winn| 
et al. ( 2010a| l, could be seen as evidence that realignment occurs 
faster than orbital decay. In a very recent development, |Lai| 



( 201 l| l shows how planets could realign their orbit faster than 
their semimajor axis would decay. 

Finding out about the ultimate fate of hot Jupiters is of 
great interest and a subject of intense on-going research, 
fraught with challenges. For example, the tidal circularisation 
and realignment timescales notably depend on the orbital 
obliquity i/r, the ratio of masses, the scaled radius (a//?*) and 
the tidal quality factors, in the planet and in the star, 
( [Hutl[T98T] [Barker & O gilvie 2009 ). Most of the theoretical 
work currently assumes constant Q values when Ogilvie & 
Lin ([2004 1 showed they depend on the tidal frequency. Lai 



( |2011| l indicates Q could vary for different physical processes. 
Similarly R^, is often assumed constant when clearly, in figure 
[Tja 1 .3 Mq star increases its radius by about 30% in about 4 Gyrs. 

Stellar age estimates are notoriously difficult to obtain. The 
estimates that have been used here have been extracted by a vari- 
ety of authors using different techniques on different sets of evo- 
lution models. The pattern resisted a blurring caused by system- 
atic effects, displaying a certain robustness. Nevertheless, this 
letter should also be an incentive to continue obtaining Rossiter- 
McLaughlin measurements as well as check those stellar ages 
and derive them in a uniform manner Similarly, accurate and 
precise age estimates for solar mass stars are dearly needed. One 
can access those via good determination of stellar parameters, 
using higher resolution spectroscopy for the T^s and Z, and high 
precision photometry which will give p*. Stellar ages can also 
be estimated from asteroseismologic timeseries underlying the 
interest in having a planet-finding space mission with such ca- 
pacity, like the proposed PLATO. Astrometric distance measure- 
ments from the GAIA satellite will soon give us an independent 
access to stellar radii. 
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Table 1. Stellar and planetary parameters used to create figure|2] Values obtained using Exoplanet.eu. The values of fi here are the 
absolute values. Error bars for the age are consistently the lower error bars presented in papers. 



Name 



PC) 



Age (Gyr) M, (Mj„p) 



references 

"(Pal et al. 2010';'Loeillet et al.'2008) 



HAT-P-2 

HAT-P-4 

HAT-P-6 

HAT-P-7 

HAT-P-8 

HAT-P-9 

HAT-P-13 

HAT-P-14 

HAT-P-16 

HAT-P-30 

HD 17156 

HD 149026 

TrES4 

WASP-3 

WASP-7 

WASP- 14 

WASP- 17 

WASP- 18 

WASP-33 

WASP-38 

XO-3 

XO-4 



1.36 
1.26 
1.29 
1.47 
1.28 
1.28 
1.22 
1.38 
1.22 
1.24 
1.28 
1.30 
1.39 
1.23 
1.28 
1.32 
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